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1. Introduction 

The analogy between quantum dots and natural atoms is rather appealing, and in 
many cases quantum dots are referred to as artificial atoms (or molecules) [1]. Within 
the physics of low-dimensional electronic systems, quantum dots and natural molecules 
realize the ultimate extreme of zero dimension. So far, most investigations have been 
focused either on planar quantum dots or on quantum dots which occupy a small volume 
(quantum box or cavity). A novel class of zero-dimensional systems which so far has not 
received much attention is realized when electrons are confined to move on a compact 
surface of nanoscopic size. The simplest class of such surfaces has the topology of the 
sphere. An electron hopping between carbon atoms of a Cgo fullerene (or its derivatives) 
provides the most natural candidate for such systems. 

In the companion work [2] we have investigated in detail the spectra of tight-binding 
electrons moving on polyhedral graphs with spherical topology, subject to the radial 
magnetic field produced by a quantized magnetic charge. This problem was solved for 
the five Platonic solids (tetrahedron, cube, octahedron, dodecahedron, icosahedron) , 
the Ceo fullerene and a couple of less symmetric objects (diamonds and prisms). The 
main goal of the present work is to pursue this idea further by including the electron 
spin and taking into account the spin-orbit interaction. Within the same framework as 
in [2], the single-particle energy spectrum of an electron subject to a radial electric field 
which generates a Rashba-type spin-orbit interaction [3] is studied for the five platonic 
solids and the C 6 o fullerene. 

The spin-orbit interaction is known to have a profound impact in atomic, nuclear 
and solid-state physics. Confining our discussion to the latter field, a dramatic example 
of its effect is the occurrence of an Anderson metal-insulator transition in disordered 
two-dimensional electronic systems [4]. Recall that the spin-orbit interaction emerges 
as a natural consequence of the Dirac equation, when the low-energy sector is described 
by the Pauli equation, and relativistic corrections are taken into account by means of 
a systematic 1/c 2 expansion [5]. For an electron of mass m and charge — e, subject to 
an electrostatic potential V(r), and therefore to a static electric field E(r) = — VV(r), 
the spin-orbit term in the effective Hamiltonian is 

eh 

Hso = t-y- 2 (p-(Ex<t) + (Ex<t)- P ), (1.1) 

am z c z 

where p = —ih'V is the momentum operator and er is the vector of Pauli matrices, so 
that S = ha jl is the electron spin operator. If the electrostatic potential V(r) = V(r) 
is central, the spin-orbit Hamiltonian (1.1) simplifies to 

f i ^_ e dV(r) L s 
2m 2 c 2 r dr 

where L is the electron orbital angular momentum operator. In particular, if the electron 
is confined to move on a spherical shell with radius R, the spin-orbit Hamiltonian (1.2) 
acquires the familiar form 

Hsq = CL-S = ^-(J 2 -L 2 -S 2 ), (1.3) 



Tight-binding spectra on spherical graphs II: the effect of spin-orbit interaction 



3 



where C is a constant and J = L + S is the total angular momentum. 

The present work will be focused on the example of the Coulomb potential produced 
by a static electric charge q placed at the center of the sphere, 

V{r) = ± E(r) = %. (1.4) 

For this potential, one has 

C = ~2rJc*R 3 - (L5) 
For the sake of completeness, we present at the end of Section 2.1 a discussion of the 
order of magnitude of the spin-orbit interaction, although this can be found in many 
textbooks. 

The Hamiltonian (1.3) has two eigenvalues E±, respectively corresponding to the 
vectors L and S being parallel and antiparallel. If £ = 0, 1, 2, . . . denotes the orbital 
quantum number, the eigenvalues E± and their multiplicities m± read 

E + = m + = 2(1 + 1), 

g _ = - C ^ + 1 >, m _=21 (1.6) 

This spectrum is not an even function of the coupling constant C, except in the classical 
regime, i.e., in the £ — > oo limit. This lack of a symmetry is expected on physical 
grounds. In the Coulomb case, C is indeed proportional to the product qe of both 
charges. Charges of the same sign (C < 0) and charges of opposite signs (C > 0) indeed 
correspond to physically distinct situations, which are not related by any symmetry. 
Furthermore, the situations where L and S are parallel and antiparallel are also known 
to exhibit different features, e.g. in scattering theory [6]. 

Our main objective is to construct and study natural discretizations of the spin-orbit 
Hamiltonian (1.3), within a tight-binding model where the electron lives on the sites 
(vertices) of a polyhedral graph drawn on the unit sphere and executes nearest-neighbor 
hopping. Our analysis will be based on an analogy with the more conventional situation 
of tight-binding (spinless) electrons subject to a given magnetic field B{r) = V x A(r). 
In this case, the hopping of particles from site A to site B is described by a hopping 
term of the form a^U abo-b + ^-c. in the tight-binding Hamiltonian, where U\b is a phase 
factor, i.e., an element of the Abelian gauge group U(l). It is generally accepted that 
the following expression, known as the Peierls substitution [7], is an appropriate choice: 

f/ AB = exp{^/ A-dr), (1.7) 

{ kc 7 7 (A,B) J 

where 7(A, B) is a given continuous path joining site A to site B. The phase factor so 
defined depends in general on the whole path 7 (A, B), and not only on the endpoints A 
and B (see [8] and [9] for recent investigations related to this matter). For the present 
problem involving the non- Abelian gauge group SU(2), the construction of the hopping 
terms requires some extra care. Indeed, since SU(2) matrices do not commute among 
themselves, an ordering prescription is needed in general. 
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The setup of the present paper is the following. The model is introduced in 
Section 2.1. The hopping matrices U AB , which are elements of the non-Abelian SU(2) 
gauge group, are evaluated for the two natural choices of shortest paths, the straight line 
segment and the arc of a great circle. A unique dimensionless parameter \x then appears 
in a natural way. The main properties of the model, and especially its symmetries, are 
studied in Section 2.2. In Section 3 the five regular polyhedra or Platonic solids and the 
C 60 fullerene (modeled as a regular truncated icosahedron) are investigated in detail. 
The spectra of the tight-binding Hamiltonian are respectively determined in Sections 3.1 
to 3.6. The total energy at half filling is studied in Section 4, whereas Section 5 contains 
a short discussion. 



2. The model 

2.1. Definitions 

In this work we consider a tight-binding model defined on polyhedral graphs drawn on 
the unit sphere. We denote by V the number of vertices (sites), by L the number of 
links (bonds) and by F the number of faces of a polyhedron. In the present case of 
spherical topology, the Euler relation reads (see e.g. [10]) 

V-L + F = 2. (2.1) 

In all the polyhedra considered in the following, all the links have equal arc length 
(with < < 7r). For any pair of neighboring vertices A and B, we thus have 

A-B = cosQ, Ax B = n AB sinB, (2.2) 

where A is the unit vector joining the center of the sphere to A, and so on, whereas 
n AB is the consistently oriented unit vector perpendicular to A and B, so as to have 

n AB = -n BA . (2.3) 

The tight-binding model is defined by means of the Hamiltonian 

H = (a{U AB a B + h.c), (2.4) 

<AB> 

where the sum runs over the L oriented links < AB > of the polyhedron, whereas 

a A = (a Ap a Al ), a A = (^^J > ( 2 - 5 ) 

where a Aa and a Aa are respectively the creation and annihilation operators of an electron 
at site A with spin component a — \ or |, and the matrices U AB are elements of the 
non-Abelian gauge group SU(2), i.e., 2x2 unitary matrices with unit determinant, 
describing the spin-orbit coupling on an electron hopping from site A to a neighboring 
site B. 

In analogy with the Abelian case described by the Peierls substitution (1.7), the 
SU(2) matrix U AB is expressed as a path-ordered integral: 



U ab = P exp -ig / (E x a) ■ dr }, (2.6) 

I y 7 (A,B) J 
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where 7(A, B) is a given path joining site A to site B, and E is the static electric field, 
as in (1.1). 

The value of the coupling constant, 

o = £?• < 2 - 7 ' 

is determined along the line of thought used in deriving the Peierls substitution in the 
Abelian case [7], and already considered e.g. in [11] in the case of the SU(2) group. 
The basic idea is to consider the spin-orbit term (1.1) as a perturbation of the free 
non-relativistic Hamiltonian H Q = p 2 / (2m) and to use the approximation 

n + n so = ^- + -—^ (p-(ex<t) + h. c .) 

~ — [P+z kE x cr) = (V + i- E x a ) . (2.8) 

2m v Amc 2 1 2m v 4m<r ' 

9 

The expressions inside parentheses in the second line are respectively the SU(2) covariant 
momentum and derivative [11]. It is worth noticing that the term obtained by expanding 
the square, i.e., g 2 h 2 E 2 / (2m), if not neglected, is a scalar potential which does not affect 
the spin physics anyhow. 

The present work is restricted to the situation where E is the electric field generated 
by a static charge q sitting in the center of the sphere, given by (1.4), so that 

(E x cr) ■ dr = -^(r x dr) • cr. (2.9) 

Let us now make the hypothesis that the path 7(A, B) is planar, i.e., entirely contained 
in the OAB plane. In this case, at every point of the path the infinitesimal vector 
r x dr is perpendicular to the latter plane, i.e., aligned with the vector ?t,ab introduced 
in (2.2). As a consequence, the path-ordering prescription is not needed, and (2.6) can 
be recast as 

%b- xp { w (/ 7(ab) ^:)..}. (2.io) 

There are two natural choices for the path 7(A, B): 
Straight-line path. If 7(A, B) is the shortest path in three-dimensional space, i.e., the 
straight line segment joining the points A and B, r can be parametrized as 

r = (l-t)A + tB (0<t<l). (2.11) 

We have then 

r x dr = sin n AB dt, r = (1 - 4 sin 2 (0/2) t(l - t)) 1/2 , (2.12) 
so that (2.10) yields 

U A b = exp(igq sin I \Q)n AB ■ cr), (2.13) 

with 

f 1 dt 1 

/(0) = Jo (l-4sin 2 (0/2)t(l-t)) 3 /2 = cos 2 (0/2)' (2 ' 14) 
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i.e., 

Uab = exp(2iggtan(6/2) n AB ■ <j). (2-15) 

Great-circle path. If 7(A, B) is the shortest path on the sphere, i.e., the arc of the great 
circle passing through A and B, r can be parametrized as 

r= sin(e-r)A + sinrB (Q < r < Q) _ (2 16) 
sin 6 

We have r = 1, whereas 

r x dr = n A B dr, (2.17) 
so that (2.10) yields 

Uab = exp(igq<~) n AB ■ a). (2.18) 

For both choices of the path 7(A, B), Equations (2.15) and (2.18) yield the same 
expression for the SU(2) matrix Uab'- 

Uab = exp(i/x Uab • tr) = cos /i + i sin/i tiab • (2-19) 

which gives the desired discretization of the familiar spin-orbit operator L ■ S recalled 
in (1.3). The fact that the hopping matrix Uab involves a vector parallel to Ax B was 
already noticed in the case of a Rashba spin-orbit interaction in semiconductors [12, 13]. 

For a given polyhedron, the model therefore has one single parameter, the 
dimensionless coupling constant /i. Re-inserting for a while the physical radius R of 
the sphere, and using the expression (2.7) of the coupling constant g, we are left with 
the following expression for //, for both choices of the path 7 (A, B): 

j2tan(0/2) (straight-line path), 
[ B (great-circle path), 

with 

qq mR 2 C qe 

e = — = = - - . 2.21 

R 2 Amc 2 R y ' 

The dimensionless (positive or negative) number e gives a measure of the strength 

of the spin-orbit interaction. Its expression (2.21) can be made more transparent by 

introducing the (positive or negative) atomic number Z, such that the charge at the 

center of the sphere is q = Ze. One has then 

£= — R> (2 - 22) 
where a = e 2 /(hc) ~ 1/137 is the fine structure constant and a = h 2 /(me 2 ) is the 
Bohr radius, whereas R is the radius of the spherical sample. Although the number e 
is a priori very small, due to the factor a 2 ~ 10~ 4 , it is allowed to become appreciable 
in the following two ways. First, there is a priori no upper limit on the value of Z, as 
the charge q = Ze is treated in this work as a static classical charge. Second, spin-orbit 
interactions can be many orders of magnitude larger in solid materials than in vacuum, 
due to Bloch electrons moving close to atomic nuclei with relativistic velocities [14]. 
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The dependence of the parameter \i on the angle 6 in (2.20) also deserves a word 
of comment. The same linear growth at small angles, i.e., fi ~ eQ, holds for both paths, 
in accord with the expectation that we are dealing with a bona fide discretization of 
the familiar spin-orbit Hamiltonian (1.3). On the contrary, the regime of large angles 
(0 — > 7r) exhibits two very different kinds of behavior: \x remains finite in this limit in 
the case of a great-circle path, whereas it diverges in the case of a straight-line path, 
as the latter passes very near the center of the sphere, where the electric field becomes 
infinitely large. 

Throughout the following, we shall adopt the theoretical viewpoint of considering /i 
as an arbitrary parameter, forgetting both about its physical origin and about its 
expression (2.20). Of course, the arc length O of the links is bound to keep its value, 
dictated by the geometry of the graph under consideration. The dependence of energy 
spectra on /i will be investigated systematically, starting with a study of its symmetries 
in the next section. 



2.2. Properties 

We now turn to a discussion of various properties of the tight-binding Hamiltonian 7i of 
the problem, defined in (2.4), where the SU(2) matrices Uab are given by (2.19), putting 
a special emphasis onto symmetries. 

Hermitian matrix representation. The relation (2.3) ensures that the matrices Uab obey 
Uab = C/ba = U BA . (2.23) 

The Hamiltonian 7i is therefore represented by a 2V x 2V Hermitian matrix 7i, whose 
rows and columns are labeled by a couple (Ac) where A = 1, . . . , V denotes a site and 
a = | or | is a spin index, such that 

Ti(Aa)(Br) = (#Ab) ctt ■ (2.2A) 

The equation for the energy eigenvalues E a , labeled by the integer a — 1, . . . , 2V, and 
the corresponding eigenf unctions ?/>A,a reads 

Ea^Ka = E U AB lP B ,a, (2.25) 
B(A) 

where B(A) runs over the neighbors of A. More explicitly, 

E a il>A*,a= E E (CabUW (2.26) 

B(A) r=U 

Sum rules. The spectrum of the Hamiltonian TL obeys the following sum rules 

E^a = 0, E S a = 4L. (2.27) 

a a 

where the sums run over the 2V eigenvalues E a , repeated according to their multiplici- 
ties. The first sum equals tr7Y = J2a^ t Uaa = 0. This sum rule is a common feature 
to all tight-binding Hamiltonians with only non-diagonal matrix elements. The second 
sum equals trH 2 = J2ab^ t (UabUba) = 4L. Equation (2.23) indeed implies that each 
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link < AB > gives two contributions equal to ^{Ux&U)^) = tr 1 = 2, i.e., the number 
of spin degrees of freedom. 

Kramers degeneracy. All the energy levels of the Hamiltonian 7i are at least 
twofold degenerate, because of time-reversal symmetry, embodied in the Kramers 
theorem [15, 16]: if ?/>A,a is a solution of (2.25), another independent solution of the 

' (K) 

same equation, with the same energy E a , is provided by the spinor V'a.oj where 

i-e., {{pi-*' (2-28) 

Here and throughout the following, the star denotes complex conjugation. 

Homogeneous modes on regular polyhedra. In the case of the five regular polyhedra, 
one can predict the existence of a twofold degenerate eigenvalue associated with 
homogeneous modes. Using the expression (2.19) of the matrices Uab, the eigenvalue 
equation (2.25) can be recast as 

Eip A = (cos pi + i sin \x n AB ■ <x)t/? B - (2.29) 

B(A) 

For each site A, consider the vector 

Wa — Ax b = sin e 2 ™ AB - ( 2 - 30 ) 

B(A) B(A) 

In the case of a regular polyhedron, one has Wa = by symmetry. Indeed, Wa is 
perpendicular to A, and in the plane perpendicular to A it has the p-fold rotational 
symmetry of the polyhedron, where p is the coordination number of the vertices. Wa 
therefore clearly vanishes. As a consequence, (2.29) shows that the homogeneous 
wavefunction t/?a = Xi where % is a constant spinor, independent of the site A, is 
an eigenfunction of the Hamiltonian Ti. The corresponding twofold degenerate energy 
is E = pcos/i. 

Semi-periodicity. The matrices Uab given in (2.19) obey Uab{^ + tt) = —Uab{^)- The 
energy eigenvalues therefore obey the same property, referred to as semi-periodicity: 
they are changed into their opposites if /x is changed to fj, + ir. It is also worth noticing 
that the spectrum of the Hamiltonian H is not an even function of /i, in spite of the 
identity Uab(—/^) = ^abIa 4 )- This lack of symmetry has already been emphasized in 
the simpler example of the spin-orbit Hamiltonian (1.3). 

//•<-> — /i symmetry. The Hamiltonian H has the following less obvious symmetry. 
For each site A, consider the spin operator in the direction of A, 

S A = A ■ a. (2.31) 
One has clearly S\ = 1. Furthermore, using the identity 

(a-tr)(b-(r)=a-b + i(axb)-(T, (2.32) 
one can check that 



Uab(Q) = S a Sb 



(2.33) 
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for any pair of neighboring vertices A and B, where ?7ab(0) is a shorthand for the 
matrix Dab given in (2.19) for \i = 0. Some algebra involving a repeated use of the 
same identity (2.32) allows one to prove the more general relation 

Uab{& -») = S A U AB (fi)S B . (2.34) 

As a consequence, the Hamiltonians H(h) and 7i(0 — //) have the same spectrum. More 
precisely, if t/?a is an eigenfunction of H(fi) with energy E, (2.34) shows that xa = Sa^a 
is an eigenfunction of H(Q — li) with the same energy E. 

Special values of li. The following two values of the parameter //: 

Ho = 6/2, /i! = 0/2 + 7T, (2.35) 

are special in several respects. The /j, «-> — /j, symmetry and the semi-periodicity imply 
that llq and Hi are symmetry axes of the energy spectrum, if displayed as a function 
of li. The Kramers degeneracy and the // <-> — // symmetry imply that all energy 
levels are at least fourfold degenerate, i.e., that all the multiplicities are multiples of 4. 

There is also a striking correspondence between the present problem at the special 
value Ho and the tight-binding problem in the presence of a magnetic monopole 
investigated in [2]. The mapping between both problems goes as follows. For li = 
Ho = 0/2 one has 

2 cos^o tW/io) = 1 + S A S B . (2.36) 

This formula, which is essentially equivalent to (2.33), suggests to introduce the local 
basis of eigenstates of the spin operators S A - Denoting by (9 a, <Pa) the spherical 
coordinates of A, such that 

A= (sin#ACOs</2A,sin#Asin<^A,cos#A), (2.37) 

the spinors Xa sucn that SaXa = ^Xa rea d 

\ / . Oa _ i¥ , A 

/ — sin p ^ A 



Xa=\ fl . I > Xa 



?A 

sm — e 



1VA 



V 



sin — e 
2 

#A 

cos ■ 



(2.38) 



2 ' v 2 

The spinors Xa are changed into one another by time-reversal symmetry, according 
to (2.28). Expanding the eigenfunctions of Tt as 

</>a = uaXa + ^aXa> ( 2 - 39 ) 

some algebra using (2.36) shows that the amplitudes ua and va obey the following two 
scalar tight-binding equations: 

Eua = J2 ^ab^b, Ev A = J2 ^b^b, (2.40) 

B(A) B(A) 

where the hopping rate £ab is given by 

2 cos^o t A B = (xt I Xb) = cos y cos y + sin y sin y e ^-^\ (2.41) 
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The expression (2.41) of £ab can be drastically simplified as follows. Using trigonometric 
identities and the relation 

cos0 = cos 9a cos 9 b + sin^Asin^B cos(<^b — V?a), (2.42) 

the expression for the square modulus of £ab can be shown to boil down to |£ab| 2 = 1- 
The hopping rates are therefore phase factors. Setting £ab = exp(iu;AB), we obtain 

cos 6 + cos 6 A + cos 6 B + 1 
cos^ab = g Ta (2.43) 

4 cos — cos — cos — 
2 2 2 

This expression can be recognized as one of the variants of the spherical Heron formula 
giving the solid angle of a spherical triangle in terms of its arc lengths [17, 18], recalled 
in the Appendix of [2]. We thus obtain 

fAB=exp(^V (2.44) 



V 2 

where ^nab is the solid angle of the oriented spherical triangle NAB, where N is the 
North pole of the unit sphere. It can be checked that phases and orientations are 
consistent, so that the product of phase factors living on the anticlockwise oriented 
links around any face equals exp(if2/2), where fl is the spherical angle of the face under 
consideration. This is precisely the requirement to describe the magnetic flux generated 
by a magnetic monopole of unit charge (n = 1) sitting at the center of the sphere. 

We have therefore shown that the spectrum of the present problem at the special 
value n — no — 0/2 consists of two independent copies of the spectrum of the magnetic 
monopole problem for n — ±1. One may wonder how the Hamiltonian Ti, which is 
invariant under time reversal, can have at the special point fj, — fj, — 0/2 the same 
spectrum as the Hamiltonian of the magnetic monopole problem, for which the time- 
reversal symmetry is broken. The key to the answer is that the expansion (2.39), which 
is the step that seemingly breaks time-reversal symmetry, amounts to performing a 
unitary transformation T> which brings the 2V x 2V Hamiltonian matrix 7i into the 
block diagonal form 

H = VH&=( H o + £ ^, (2.45) 

where 7i± are the V xV Hamiltonian matrices of the magnetic monopole problem with 
respective magnetic charges n — ±1. The full Hamiltonian 7i is invariant under time 
reversal. The operators 7i + and = 7i+ separately break time-reversal symmetry, 
but they are changed into one another by time reversal. Each of them brings one copy 
of the spectrum of the magnetic monopole problem. 



3. Polyhedra and their spectra 

In this section we investigate the spectrum of the tight-binding Hamiltonian 7i for the 
five regular polyhedra or Platonic solids and for the fullerene, modeled as a symmetric 
truncated icosahedron. In the case of the Platonic solids, all the properties derived in 
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Section 2.2 will be checked against our analytic expressions of the energy spectra. The 
fullerene will be special in the following two respects. Its energy spectrum will not be 
obtained analytically, albeit from the numerical diagonalization of an explicit 120 x 120 
matrix. The homogeneous modes described in Section 2.2 are absent, as the fullerene is 
not sufficiently symmetric to allow them. 

All these polyhedra have been described in detail in [2]. Table 1 lists a few of their 
geometrical characteristics which will be useful in the following. 



polyhedron 
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L 
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tetrahedron 
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-1/3 


cube 


8 


12 


6 


3 


4 


1/3 


octahedron 


6 


12 


8 


4 


3 





dodecahedron 


20 


30 


12 


3 


5 




icosahedron 


12 


30 


20 


5 


3 




fullerene 


60 


90 


32 


3 


C 


(80 + 9^)7109 



Table 1. Geometrical characteristics of the polyhedra considered in this work: 
numbers V of vertices, L of links, F of faces, coordination number (number of neighbors 
of a vertex) p, number of sides of a face q, expression of cos0, where the arc length 6 
of the links has been introduced in (2.2). 

3.1. The tetrahedron 

The tetrahedron is the simplest of the Platonic solids. It consists of 4 trivalent vertices, 6 
links and 4 triangular faces. 

Throughout the following it will be advantageous to unwrap the polyhedra around 
an axis of high symmetry, say of order r, to be used as the z-axis [19]. For all the Platonic 
solids, the order r of rotational symmetry can be chosen to be the larger of the integers p 
and q. The planar representation of the tetrahedron thus obtained, emphasizing the 
vertices and the links between them, is shown in Figure 1. Some vertices and links may 
have several occurrences, to be identified by the inverse procedure of wrapping the planar 
representation onto the sphere. This planar representation is an efficient tool to find the 
Cartesian coordinates of the vertices, making an optimal use of symmetries. Vertices 
at the same height on the plot have the same z coordinate, whereas their coordinates 
in the xy-plane are obtained from each other by rotations by the commensurate angles 
lukjr for k = 1, . . . , r. Table 2 lists the coordinates of the vertices thus obtained. 

The 8x8 Hamiltonian matrix 7i has been constructed from these coordinates by 
using (2.19) and (2.24). This matrix is too complex to be diagonalizable by hand. 
For each of the five Platonic solids this task has been performed with the help of 
the software MACSYMA. The energy eigenvalues and their multiplicities are listed 
in Table 3. Horizontal lines separate groups of levels related to each other by the 
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Figure 1. Planar representation of the tetrahedron. 



vertex 


X 


y 


z 


A 








1 


B 


2V2/3 





-1/3 


C 


-V2/3 




-1/3 


D 


-V2/3 


-V6/3 


-1/3 



Table 2. Cartesian coordinates of the vertices of the tetrahedron. Horizontal lines 
separate groups of vertices having the same z coordinate. 

//<->© — // symmetry. The levels Ei(ff) and E 2 ((f) are interchanged by this symmetry, 
whereas E 3 (/i) is symmetric by itself. One has indeed 

E 2 { f i)=E 1 {Q-fi), EM = -V3cos(ii-e/2). (3.1) 

The energy spectrum is shown in Figure 2 as a function of \ij (2ir) over one period. The 
vertical dashed lines show the symmetry axes of the spectrum at the special values of \i 
given in (2.35). 



a 


£a(/i) 


m a 


1 


3 cos/i 


2 


2 


— cos /i + 2v^2 sin jx 


2 


3 


— cos ji — \/2 sin \x 


4 



Table 3. Energy levels E a (p) of the tetrahedron and their multiplicities m a . 
Horizontal lines separate groups of levels related to each other by the \x <-» 6 — yU 
symmetry. 

3.2. The cube 

The planar representation of the cube is shown in Figure 3. Table 4 lists the Cartesian 
coordinates of the vertices. 

The energy eigenvalues of the 16 x 16 Hamiltonian matrix constructed from these 
coordinates are listed in Table 5 and shown in Figure 4 as a function of ix/{2n) over 
one period. The spectrum is observed to be its own opposite, i.e., to be symmetric with 
respect to the origin of energies, E — 0. This extra symmetry is particular to the cube, 
being due to the fact that this polyhedron is bipartite. 
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3 
2 

1 

H 
-1 

-2 

~ 3 0.2 0.4 0.6 0.8 1.0 
n/(2n) 

Figure 2. Plot of the energy spectrum of the tetrahedron against fi/(2n) over one 
period. Vertical dashed lines: symmetry axes at the special values of \x given in (2.35). 

A B C D A 
T T T 1 T 

i i 1 i i 

E F G H E 
Figure 3. Planar representation of the cube. 
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-i/Vs 


-i/Vs 
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i/Vz 
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-1/V3 
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-i/Vs 
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-l/v/3 


H 




-1/V3 


-l/v/3 



Table 4. Cartesian coordinates of the vertices of the cube. Same conventions as in 
Table 2. 

3.3. The octahedron 

The planar representation of the octahedron is shown in Figure 5. Table 6 lists the 
Cartesian coordinates of the vertices. 

The energy eigenvalues of the 12 x 12 Hamiltonian matrix constructed from these 
coordinates are listed in Table 7 and shown in Figure 6 as a function of fx/ (2ir) over one 
period. 
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/ > i 
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~.j y y j i j lis 
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cos fx + 2v^2 sin fx 


2 
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cos fx — y/2 sin fx 
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— cos /i + sin /x 


4 


5 


— cos fx — 2\/2 sin fx 


2 


6 


—3 cos /i 


2 



Table 5. Energy levels E a {ix) of the cube and their multiplicities m a . Same 
conventions as in Table 3. 




0.2 0.4 0.6 0.8 1.0 



Figure 4. Plot of the energy spectrum of the cube against fx/(2w) over one period. 
Same conventions as in Figure 2. 



A A A A 




F F F F 
Figure 5. Planar representation of the octahedron. 



3.4- The dodecahedron 

The planar representation of the dodecahedron is shown in Figure 7. Table 8 lists the 
Cartesian coordinates of the vertices, with the shorthand notations c& = cos(/c7r/5), 
s k = sin(A;7r/5), and 
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Table 6. Cartesian coordinates of the vertices of the octahedron. Same conventions 
as in Table 2. 



a 




m a 


1 


4cos/i 


2 


2 


4 sin \i 


2 


3 


—2 sin fi 


4 


4 


—2 cos \x 


4 



Table 7. Energy levels E a (p) of the octahedron and their multiplicities m a . Same 
conventions as in Table 3. 




0.2 0.4 0.6 0.8 1.0 
/V(2tt) 



Figure 6. Plot of the energy spectrum of the octahedron against /u/(27r) over one 
period. Same conventions as in Figure 2. 




Figure 7. Planar representation of the dodecahedron. 



The energy eigenvalues of the 40 x 40 Hamiltonian matrix constructed from these 
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Table 8. Cartesian coordinates of the vertices of the dodecahedron. Same conventions 
as in Table 2. Shorthand notations are explained in and above (3.2). 



coordinates are listed in Table 9 and shown in Figure 8 as a function of /i/(2rr) over 
one period. This is the first example where some of the energy levels are not given 
by linear functions of cos/x and sin/x. The four sixfold degenerate energy levels £j(/x) 
(i = 1, 2, 3, 4) are the roots of the polynomial equation 

£ 4 + 2£e 3 + (2? - 5)e 2 + 2£(£ 2 - 4)e + 1 - 3£ 2 + £ 4 = 0, (3.3) 

with 

£ = cos /x + sin/x = \/3 cos(/i — 0/2), (3.4) 

and where the branches are chosen such that £i(0) = 1, £2(0) = 0, £3(0) = —2, 
£ 4 (0) = -y/b. 



a 
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3 cos ji 


2 
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y/h cos /i + 2 sin /x 
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3 


\/5 cos /x — sin /x 


4 


4 


cos /x + y/b sin /x 
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5 


£!(/x) 


6 


6 


£ 2 (/i) 


6 


7 


—3 sin fi 


2 


8 


—2 cos /x + \/5 sin /x 


2 


9 




6 


10 


£ 4 (/i) 


6 



Table 9. Energy levels E a (fi) of the dodecahedron and their multiplicities m a . Same 
conventions as in Table 3. The functions £j(/x) (i = 1,2,3,4) are the roots of the 
polynomial equation (3.3). 
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/V(2tt) 

Figure 8. Plot of the energy spectrum of the dodecahedron against ^/(2tt) over one 
period. Same conventions as in Figure 2. Thick dotted lines: sixfold degenerate energy 
levels £j(/u) (i = 1,2,3,4). 



3.5. The icosahedron 

The planar representation of the icosahedron is shown in Figure 9. Table 10 lists the 
Cartesian coordinates of the vertices, with the shorthand notations Ck = cos(kn/5), 
Sfc = sin(A;7r/5), and 

d=M, * = (3.5) 



A A A A A 




L L L L L 
Figure 9. Planar representation of the icosahedron. 

The energy eigenvalues of the 24 x 24 Hamiltonian matrix constructed from these 
coordinates are listed in Table 11 and shown in Figure 10 as a function of /i/(2ir) over 
one period. 

3. 6. The fullerene 

We now turn to the case of the Cqq fullerene. For simplicity we model this molecule as 
a symmetric truncated icosahedron, where all the links have equal lengths, so that the 
analysis of Section 2.1 applies. Let us however recall that this symmetry is known to be 
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Table 10. Cartesian coordinates of the vertices of the icosahedron. Same conventions 
as in Table 2. Shorthand notations are explained in and above (3.5). 
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5 cos/i 


2 


2 


y^cos /i + 2 sin jx) 


2 


3 


V^cos fi — sin fx) 


4 


4 


— cos fx + 3 sin fx 


4 


5 


— cos fx — 2 sin fx 


6 


6 


— v^5 cos /i 


6 



Table 11. Energy levels E a (fi) of the icosahedron and their multiplicities m a . Same 
conventions as in Table 3. 




0.2 0.4 0.6 0.8 1.0 
/V(2tt) 



Figure 10. Plot of the energy spectrum of the icosahedron against fx/(2w) over one 
period. Same conventions as in Figure 2. 

slightly violated [20], as for the free molecule the length of the sides of the pentagons is 
1.46 A, whereas the length of the other links is 1.40 A. 

The symmetric truncated icosahedron has V = 60 equivalent vertices, L = 90 
equivalent links, and F = 32 faces, namely 12 pentagons and 20 hexagons, respectively 
corresponding to the vertices and to the faces of the icosahedron. Figure 11 shows the 
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planar representation obtained by unwrapping the fullerene around a fivefold axis going 
through the opposite pentagonal faces A 1 . . . A 5 and Li . . . L 5 . 




Figure 11. Planar representation of the fullerene (symmetric truncated icosahedron) . 

The Cartesian coordinates of the vertices of the fullerene have been derived 
from those of the vertices of the icosahedron, listed in Table 10, using the approach 
described in the Appendix of [2]. This procedure is illustrated in Figure 12, showing an 
enlargement of the upper left part of Figures 9 and 11, with consistent notations. One 
has 

Ai = \{2A + B), A 2 = X(2A + C), (3.6) 
and so on, with 



so that 




B B 2 C 5 C 



Figure 12. The triangular face ABC of the icosahedron decorated by vertices of the 
fullerene. Notations are consistent with Figures 9 and 11. 

The energy eigenvalues of the 120 x 120 Hamiltonian matrix constructed from the 
coordinates thus obtained have been evaluated by means of a numerical diagonalization. 
The energy spectrum is shown in Figure 13 as a function of n/(2ir) over one period. 
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For /J, — 0, i.e., in the absence of spin-orbit coupling, we recover two independent 
copies of the known tight-binding spectrum of the fullerene [21], with its 15 distinct 
energy levels with multiplicities ranging from 1 to 9. For generic non-zero values of fi, 
the spectrum consists of 28 distinct energy levels with multiplicities ranging from 2 to 6 
only. As fi — * the 28 levels merge into the 15 ones according to the patterns given in 
Table 12. We have introduced the shorthand notation 

w± = ^2(19 ± Vh). (3.9) 




\±l (2tt) 

Figure 13. Plot of the energy spectrum of the fullerene against fi/(2n) over one 
period. Same conventions as in Figure 2. 



4. Total energy 

An interesting illustration of the above energy spectra is provided by the total energy 
at half filling, defined as 

£ = X>a, (4.1) 
a=l 

where the 2V energy levels are assumed to be in increasing order (Ei < E 2 < . . . < E 2 v) 
and repeated according to their multiplicities. 

The first of the sum rules (2.27) implies that the total energy thus defined 
is insensitive to the sign of the Hamiltonian 7i. Combining this feature with the 
symmetries derived in Section 2.2, we conclude that S(/j.) obeys the symmetries 

S{n) = S{Q-n) =£(fi + n). (4.2) 

The total energy therefore has period tt, and exhibits two inequivalent stationary points 
per period, at 

^0 = 6/2, /i m = (0 + tt)/2. (4.3) 
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Table 12. Energy levels £'(0) of the fullerene at fi = 0, with their multiplicities m(0) 
and degeneracy-lifting patterns at small \i ^ 0. The shorthand notation w± has been 
introduced in (3.9). 



The first of these values, fiQ, coincides with one of the special values introduced in (2.35), 
i.e., one of the symmetry axes of the spectrum. The second of the above values, /i m , 
corresponds to one of the midpoints between the latter symmetry axes. 

The second of the sum rules (2.27) implies that the mean squared value of the 
individual energy levels is (E 2 ) = AL/(2V) = p, where p is the coordination number of 
the vertices. This suggests to introduce the reduced total energy 

- £ 8 8 (4.4) 



2V s fp 2V2VL' 



This heuristic argument can be turned to a quantitative prediction in the p — > 00 limit 
of a very highly connected structure [2]. In this limit, the reduced total energy £ r has 
been shown to have the universal limiting value 

£co = ]= = -0.398942. (4.5) 

v 2n 



Figure 14 shows a plot of the reduced total energy £ r for all the polyhedra 
investigated in this work, as a function of (/i — yU ) / vr over one period. The reduced 
total energy is observed to wander around the limiting value (4.5), shown as a dashed 
line. The amplitude of the oscillations, i.e., of the dependence of the total energy on 
the parameter /i, is a decreasing function of the number of vertices. Figure 15 shows an 
enlargement of the plot focusing on the weak /z-dependence of £ T in the two examples 
with the larger numbers of vertices, i.e., the dodecahedron (V = 20) and the fullerene 
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-0.25 




' 0.2 0.4 0.6 0.8 1.0 
(/x-/i )/rr 

Figure 14. Plot of the reduced total energy £ r against (fi— ho)/tt 1 for all the polyhedra 
investigated in this work: tetrahedron (V = 4) (empty diamonds), cube (V — 8) (down 
triangles), octahedron (V = 6) (up triangles), dodecahedron (V = 20) (squares), 
icosahedron (V — 12) (circles) and fullerene (V = 60) (stars). The curves consist of 
many more points than symbols (500 data points for each polyhedron) . The horizontal 
dashed line shows the limiting value (4.5). 




Figure 15. Enlargement of Figure 14, with the same conventions, emphasizing the 
weak /i-dependence of the reduced total energy E x in the cases of the dodecahedron 
and of the fullerene. 



{V = 60). The abscissa axes in Figures 14 and 15 are such that the stationary point /j,q 
introduced in (4.3) corresponds to the ends of the plots, whereas fi m corresponds to 
their centers. The latter stationary point is observed to be the absolute minimum of 
the total energy for all the polyhedra considered in this work, except the tetrahedron, 
for which the total energy has its absolute minimum at // = /i and a local minimum 
at n = /i m . 
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5. Discussion 

In this paper we have introduced and investigated a tight-binding model defined on 
graphs drawn on the unit sphere, describing the motion of an electron subject to a 
spin-orbit interaction in the radial electric field created by a classical charge sitting at 
the center of the sphere. The present work completes our study of electronic properties 
of mesoscopic and nanoscopic systems with the topology of the sphere, started in the 
companion work [2], which is devoted to electrons subject to a radial magnetic field 
produced by a quantized magnetic charge sitting at the center of the sphere. 

This work has been focused onto polyhedral graphs such that all links have a 
common arc length 0. For a fixed graph of this kind, the model has only one 
parameter, /i, giving a dimensionless measure of the strength of the spin-orbit inter- 
action. Among the symmetry properties of the model, exposed in detail in Section 2.2, 
the /i <-> — /i symmetry was quite unexpected, as it has no counterpart in the 
continuum, described by the familiar form L ■ S of the spin-orbit Hamiltonian, whose 
eigenvalues and multiplicities are recalled in (1.6). 

For the special value /i — /io — 0/2, which coincides with one of the symmetry axes 
of the spectra, an exact correspondence has been established with the tight-binding 
problem in the magnetic field of a Dirac monopole, investigated in [2]. It is remarkable 
that, by tuning a parameter in a theory which is experimentally realizable, it is possible 
to obtain the spectrum of another system whose experimental realization is so far elusive. 
In fact this correspondence already holds at the classical level. It was indeed discovered 
long ago by Poincare [22] that the motion of an electrically charged particle in the field of 
a magnetic charge can be mapped onto that of a spherical top. In a quantum-mechanical 
framework, the quantitative correspondence between both problems reads [23, 24] 

1711 = 25, (5.1) 

where the integer n is the magnetic charge (in units of the elementary magnetic charge 
of Dirac's monopole), whereas S is the total spin of the top. The above relation can 
easily be recovered by noticing that the ground state of the Schrodinger equation on 
the sphere in the presence of a magnetic charge n, investigated in the pioneering work 
of Tamm [25], has a multiplicity \n\ + 1, to be identified with 2S + 1. In the present 
situation of an electron (S = 1/2), the correspondence has indeed been shown to hold 
for a unit magnetic charge (n = ±1). More generally, a similar correspondence can 
be expected to hold true for higher representations as well, whenever the spin and the 
magnetic charge are related by (5.1), for a suitably chosen special form of the spin-orbit 
interaction. 

We have then turned to the study of specific examples of polyhedra, namely the five 
Platonic solids (tetrahedron, cube, octahedron, dodecahedron and icosahedron) and the 
C 60 fullerene (modeled as a regular truncated icosahedron). For the Platonic solids, the 
full /^-dependence of the energy levels, and the corresponding multiplicities, have been 
obtained analytically in Sections 3.1-3.5. These results allow for an explicit check of the 
general properties listed in Section 2.2. Rather surprisingly, all the energy levels can be 
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expressed as linear combinations of sin// and cos/x, except the four sixfold degenerate 
levels of the dodecahedron, which are obtained as the roots of a fourth-degree polynomial 
given in (3.3). This simplicity of the energy eigenvalues is to be contrasted with the 
rather large dimension, 2V, of the Hamiltonian matrices, i.e., 40 for the dodecahedron. 

Pursuing along the lines of our companion work [2] , we have also evaluated the total 
electronic energy S of the system at half filling. For all the examples considered in this 
work, this total energy is found to be rather close to its asymptotic value in the limit 
of large coordination numbers, where the density of states becomes Gaussian. Finally, 
as far as its dependence on the parameter /x is concerned, the total energy reaches its 
absolute minima at the midpoints between the symmetry axes of the spectra for all the 
polyhedra considered in this work, except for the tetrahedron where the total energy 
has its absolute minima at the symmetry axes of the spectrum and local minima at the 
midpoints. 
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